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We present a study of the motion that a massive particle in a non-dissipative Brownian motion 
in vacuum is subject. Noise source comes from vacuum fluctuations of some quantum field capable 
of interact with the particle. For the associated Fokker-Planck equation we do a perturbation 
theory which has terms presenting dynamics similar to that satisfied by the probability amplitude 
in Quantum Theory even though the physical interpretation of these terms be classical. In particular 
observables present expected values that coincide with those calculated using Quantum Mechanics. 
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A. Introduction 



The probabilistic character of Quantum Mechanics as well as the formal similarity of the Schroedinger equation to 
the diffusion equation raised from long ago the question of a possible stochastic justification to Quantum Mechanics 
. A first formal approach was developed by Nelson ||^ and up to now many authors presented complementary or 
alternative works (a few of them listed in [^) on the subject named usually as Stochastic Mechanics. A common 
goal of these theories lie on the possibility of find a stochastic process that is equivalent to Quantum Mechanics. 
Nevertheless no such a process was yet proposed that lent to Stochastic Mechanics, in any of its presented forms, a 
explicit advantage over the sophisticated epistemological structure of Quantum Theory. On another way Quantum 
Theory is continuously tested at deeper levels fixing its place as fundamental one although this do not invalidate the 
search for alternative interpretations, hopefully not subject to inconsistencies (like infinities) yet plaguing the theory. 
In this work we focus on a somewhat different approach where a particular kind of stochastic process is studied. 
We present a perturbation theory for the associated probability density function that has quantum formalism as a 
subjacent dynamics but still preserving the classical nature of the theory. We also discuss a possible physical origin for 
the noise source capable of describe the type of process we study. No full equivalence to Quantum Theory was found 
so this indicates that Quantum Mechanics as a fundamental theory is preserved also. Nevertheless our formulation 
shows explicitly that quantal aspects may appear quite naturally in a classical theory. 

Once every physical process include or can generate some sort of noise the motion of a massive particle under its 
influence is an important subject in the study of real particle's dynamics. For a particular class of noise - the Wiener 
process - a rigorous probability formalism can be obtained where a field equation - the Fokker-Planck equation (FPE) 
- for the probability density function of some noise-dependent variable is defined 0. For this case noise source has 
an uncorrelated stationary probability distribution and when the stochastic process under study describes specifically 
the motion of a particle the resulting process is the Brownian motion. Due to its generality a FPE can be considered 
as a fundamental although its detailed predictions are dependent of the specific Wiener process considered. 

The origin of noise that a particle in Brownian motion is subject usually is thermodynamic. This is the case of 
a particle dispersed on a fluid at some temperature. There thermal noise gives a bit of energy to the particle and 
dissipative effects returns it back (randomly) to the thermal bath. The net effect is a random walk equilibrated to the 
bath's temperature. Here we will consider this kind of motion but with two main differences. First the particle is in 
vacuum and noise is generated by vacuum ffuctuations of some quantum field capable of interaction with that particle. 
Second the motion is not dissipative. This last hypotesis seems to be at first sight unrealistic because particle may 
gain energy from nothing. We accept it for simplicity since this hypotesis do not invalidate the methods we present 
here for the general case. 

One may argue against the reality of vacuum fiuctuations itself. This subject has yet some controversy but its 
theoretical as well as experimental grounds has been sufficiently discussed in the literature |^ and we can accept 
it in a general basis. Moreover recent careful measurements of Casimir forces demonstrated a good agreement 
between observation and theoretical predictions reinforcing our belief in its physical existence. Generally we even can 
believe that geometry changes in a physically confined systems may induce a small violation of energy conservation 
principle, undetectable by today's instrument sensitivity or be possible that during molecular conformational changes 
some energy comes from vacuum, involving van der Walls force. We think these arguments are sufficient to justify 
a possible existence of a natural noise source interacting to any elementary particle coming from vacuum properties. 
Besides all we must also have in mind that every charged particle is subject to the cosmic background radiation 
field that is real and also stochastic. Therefore we think that stochastic process are inevitable to the motion of 
massive particles in vacuum, coming as a intrinsic part of the reality in Nature and must be considered in any precise 
description of this kind of motion. 

The method we give here to solve the FPE for the Brownian motion has general validity and no fundamental 
difficulties are present when the more general case where not only dissipative effects as well as a larger class of 
Wiener process are considered. Specifically for a quantum noise coming from vacuum fluctuations we understand a 
stationary Wiener- type source with noise intensity proportional to a "vacuum power " which we write asV = jh {yj^) 
and with mean energy £ = ^H^/Juj^ and or generally 

Vh^S^ (1) 

Here H is the Planck's constant and (w^) is the variance of the field frequencies averaged over some appropriate 
distribution (we assume (lo) = since uj and —u must be considered as independent fiuctuations). For example in 
the cosmic background case where T ~ 2K we find V ~ l.lbpW. Calculation of (w^) for quantum ffuctuations is 
not trivial because vacuum energy density diverges as oj^ ^ with (assumed) uniform probability distribution denying 
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a simple averaging process unless a physical cutoff at high frequencies exist. This must be achieved by multiplying 
the instantaneous vacuum modes by the appropriate cross section (in that frequency) which falls off properly at high 
energies. If doing so the instantaneous momentum transferred to the particle may be calculated resulting in a noise 
source with zero mean and the desired power in as much as in the Casimir effect calculation where high frequencies 
modes cancel out and the net effect comes from low frequencies modes only. This is analogous to the thermal effect 
where the enormous kinetic energy of one mol of molecules at room temperature has minimal action over the center of 
mass motion of a macroscopic body because the random collisions between the molecules and the body are isotropic 
so the net resulting effect are small fluctuations. Nevertheless if the phenomena we describe here has to be universal 
it is reasonable that £ should be a characteristic mensurable constant in Nature even though we concentrate here only 
on the effect on particle motion of a noise source with properties listed above and explicit in the perturbative series 
expansion procedure we give below. Therefore besides eqn detailed physical origin of V (or £) is not relevant to 
our calculations in the present stage of development of our reasoning. 



B. Perturbation expansion for the Fokker-Planck equation 



In sequence we give the precise mathematical formulation of the kind of motion we want to describe. Particle's 
dynamics arise from Hamilton equations plus a Langevin term acting as a stochastic force. We assume no explicit 
stochastic term for position variable meaning that any random effect on it comes indirectly through coupling with 
the momentum term resulting in a non-dissipative Brownian motion. More precisely we have 



ax (t) = -^dt 
op 

dH 

dp (i) = — 7^dt + \/2mVdW (t) 
Ox 



(3) 
(4) 



Here m is particle's mass and dW (t) a stationary Wiener-type stochastic variable with noise power V given as above. 
Rigorously speaking if the particle has a electric charge it radiates when accelerated which is a kind of an intrinsic 
unremovable dissipative effect. In this work we assume these effects as small thus leaving out any radiation reaction 
terms and deserving for a future work a proper inclusion of appropriate corrections. 

The associated Fokker-Planck equation for the joint $(x,p, t) probability density distribution in phase space is 
given by 



{x,p,t) 
dt 



dH a$ 

dp dx 



dHd<i> ^^92$ 
dx dp dp^ 



(5) 



which unless a missing dissipative term has the same structure of the Kramers' equation orginaly written to 
describe the motion of a large mass particle in a fluid at thermal equilibrium with a bath in such a way that noise 
comes from the temperature effect cited above. As already discussed this motion is unstable against infinitesimal 
stochastic perturbations. In fact even if the intensity y/lmV of the Langevin term goes to zero a smeared phase space 
still results. We observe that in this case the above equation reduces to the Liouville equation which admits a class 
of solutions concentrated along the deterministic trajectories in phase space describing particle's classical motion. In 

fact if H^i {t) = "l^^ -I- V {xci (t)) is the Hamiltonian of the (deterministic) Hamilton equations one class of solutions 
for Liouville equation is 



(f> {x,p,t) = Aexp(-/3 



P^ - Pel {tf 



2m 



Vix)-Vixci (t)) 



(6) 



We see that trajectories (in phase space) no longer exist any more: position and momentum variables decouple and 
information about the system is supplied by an ensemble of states in thermal equilibrium at temperature f3~^. This 
temperature effect is the only memory of the ghost stochastic effect and is caused by the known singular small noise 
expansion of the FPE §. Thus there is a big difference if we take Hamilton equations without noise to get a full 
deterministic motion and a FPE with null noise term. In the first case the topology of phase space as represented by 
Poincare maps are sets of points even when the motion is chaotic whereas a probability density function continuously 
fills phase space in the another case. 
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The structure of the FPE allows easy calculation of some averages. For example time derivative of the expected 
value of the Hamiltonian is given by 



^{H)= I dxdp 



dt 



d ^ ,dV p , ,dV p , „ 

{H)=V-{—^ + { — ^)^r (7) 
dt ox m ox m 

and the expected value of the momentum changes at the rate 

following in mean the second Newton's law. Note that the mean exchanged momentum with stochastic source is 
zero but a quadratic dispersion on momentum exists being this consequence of the missing dissipative effects. This 
dispersion increases with time meaning that, as anticipated above, the particle gains some energy from the stochastic 
source . Position variable is also affected although there is no noise source acting directly on it. Its mean value and 
the associated fluctuations change at the rate 

= (-> (9) 
dt m 

±({:c^)-{xr) = l({px)-{p){x)) (10) 
dt ' ni 

Thereby classical laws of motion are satisfied only in mean no matter the value of the noise intensity as shown by 
eqns (||)and (9). We must conclude that a field dynamics for the probability density should exist even for zero noise 
power and it will be shown bellow it follows exactly quantum dynamics when our perturbative schema is considered. 

We develop in sequence the perturbative expansion. To do this we shall use h in order to Fourier Transform 
$ {x,p,t) in momentum to the adjoint space y associated to this variable. In the present formalism this space refer 
to wavelengths of a mode decomposition of the probability density function. Differently from Quantum Theory here 
position and momentum are not related by an unitary transformation between adjoint spaces. This happens for 
momentum and wavelengths of waves of probability density having momentum p = hy. It is in this space the FPE 
has a form suitable to the structure of the perturbative series we look for. Nextly we write 

$(x,p,t) = — / x{x,y,t)e''^dy (11) 



dtxix,y,t) ^ —dyd,,x- ——X- —X (12) 



so the FPE reads 

im in ox TA'^ 

where A = £/hc is the mean vacuum fluctuation wavelength and r the particle's Compton length divided by the 
velocity of light. We stress that phase space probability density $ {x,p,t) is similar but has to be not confused with 
the Wigner function in Quantum Theory. The former satisfies exactly the Liouville equation for zero noise power. 
Wigner function satisfies it only approximately. Moreover in a Fourier transform on momentum variable of the Wigner 
function similar to eqn(|lT|) variable y has the meaning of position. Specifically in the Wigner function the role of x 
and y are symmetric so they share the same physical meaning whereas in this work simbols x and y are truUy indepent 
variables that play different role. It will be shown explicitly later that in fact both functions do not coincide. 

Owing to the structure of the eqn(|l^) it is natural to think in a Taylor expansion of x i^j Vi t) in the variable y 
which should be convenient because this way we may compare probability wavelengths to vacuum wavelength (that 
is very big) as expressed by the fraction present in the last term of this equation. This perturbative schema has the 
advantage that it permits a controU of noise effects over field dynamics for x- Saying differently if particle's kinetic 
energy is much greater than vacuum energy £ its typical mode wavelength y ~ h/p is very small compared to vacuum 
wavelength A thus making noise effects negligible. But a infinite series in powers of y has no simple inverse Fourier 
representation generating a somewhat pathological phase space reconstruction. For vanishingly small noise we know 
from eqn(^ that Liouville equation feels the presence of a thermal bath with temperature In the present case 

this effect is to be replaced by a "vacuum temperature" proportional to A^^. Thus we insert a additional gaussian 
term in the perturbation series in order to regularize phase space reconstrution and to provide the appropriate thermal 
bath effect. More precisely we look for a general series solution of the type 



2 •■ 



X(x,2;,i)=exp( ^ 1 ^a„(x,i)y" (13) 

n— 
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in such a way that eqn(|lj) becomes 

^ dUn {x, t) 



dt 



1 m ' ^ 



dan {x, t) 
dx 



ny 



..n-i y 



n+l 



A2 



ih dx 

generating the following set of differential equations for the coefficients {a„} 

dao (x, t) h dai (x, t) 



dt 

dai {x, t) h 
dt im 



im dx 

da2 {x, t) 1 daQ [x, t) 



dx 



A2 dx 



}_dV_ 

ih dx 



a-o {x,t) 



(14) 
(15) 



dttn jx, t) 

dt 



(n+l) 



dan+i (x,t) 1 9a„_i (x,t) 



9x 



A2 



dx 



(16) 



1 9y , , 1 , , 

-— «„_i(x,t)-;^a„_,(.,i) 



Before discussing solutions of the above set of equations note that a naive interpretation can be given to some of 
these coefficients. In fact since 



/oc 
$(a;,p,t) e-'^^dp 
-oo 



we have 



ao {x,t) =x{x,0,t) 
oi (x,t) : 



$(a;,p,t)dp= |*(x,t)r 



9x\ 1 f°° 1 

— = -z / p<i>{x,p,t)dp= —Tr{x,t) 



y=0 



(17) 



(18) 



that is ao is the marginal probability distribution, ai is proportional to the expected value of the momentum at point 
X and 02 is proportional to the kinetic energy at x. We have also introduced the probability amplitude \1/ belonging 
to the stuff of the L^-based probability distribution oq. Integrating eqns(14)-(16) we get for n = 2 



jJj^MUx^-=, 

d r , N , 

— / TT (X, t) dx 

dt J^oo 

-^1 fa. (.,.)..= ' <^{p')^ifM^i^^^)) 
m dt J^^ 



2m dt 



V 



The first of these equations indicates that the norm of 5* is time-independent. Incidentally it also shows that the 
Hilbert space operator satisfy {ihdt^,"^) = {'^jihdt'^) where (•, •) stands for Hilbert space inner product. This 
does not mean at the present stage that ihdt is hermitian since the relation was not proved to be valid in general, for 
any pair of Hilbert space vectors. However it is a necessary condition suggesting that this operator admits a Hermitian 
representation on an appropriate Hilbert space and in fact it will be constructed as the consistence of our calcualtions 
becomes closed. The second equation is similar to the Ehrenfest theorem and the third one indicates together with 
eqn(7) that the rise in total particle's energy comes from its kinetic component. It's worthwhile note that the first 
two equations are independent of noise intensity so concerning the non-stochastic limit their validity has the same 
status as the Liouville equation. In the next step we look for a field equation for 5*. Define L = F[P]-\- y (x,t), 
where P = j-^ .We then have 
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i ox 



and in consequence 



^ { t) ^ ^ 
dt ' im dx 



2a2 {x,t) - — ao {x,t) 



** [L, P] * 
fi2 



(19) 



d f°° 1 

-j a^{x,t)dx=-(^,[L,P]^) (20) 
Being ai pure imaginary we can write it as 

«i = ^ (21) 

where B is necessarily a real-valued operator. Further development on defining a more specific form for this 
operator also generates the rules for the Hilbert space we are working on. Owing to eqn(pO|) where a inner product 
emerges naturally we choose B such that B = iB is Hermitian and time- independent giving the desired structure to 
expected value of the momentum as shown below 

J ai {x, t)dx = ^ (^*, B*) ^ B^-) = (tt) EE n (t) 

This result is to be interpreted classically even though it presents a strong resemblance to quantum mechanical rules. 
The average force on the particle is then equal to 

dn{t) d f°° 1 



dt 



ai{x,t)dx^^ {ihdt-^ , B-^^ - , Bindt-^^ (22) 



and when compared with eqn(^0|) suggest a simple inner product structure for the above equation giving to ihdt a 
possible Hermitian representation TL which in general is function of x, P, B and t. This option defines simultaneously 
a field equation for the probability amplitude since now we select those Hilbert space vectors that satisfy the equality 

ihdt'H = n (x,P,B,tj * 

When this equation is inserted in eqn(^2|) and compared to eqn(^o|) we get 

[h,^ *) - (*,[L,P]*) 

which should be valid for all 4* that satisfy the above field equation. So in an appropriate subspace we must have 



n.B 



[L, P] with a formal solution given by 



B = P 

Ti:* = = [F [p] + V {x, t)] * 

which must be simultaneously valid. This way we conclude that F [P] is Hermitian. Now eqn(14) should also be 
satisfied implying that 

or 

(f [p] + ^di\ * = (**)-i (f [p] + ^di\ 
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an equation locally valid. Since 4" and are independent we have F [P] = ^-^d"^ and Ti. = 



(23) 



We have found that field equation for the probability amplitude \1/ (x, t) follows quantum dynamics although the 
formalism used in this work be classical. This classical signature of Stochastic Mechanics as presented here is not only 
conceptual but has rigorous formal support as shown below. 

The whole perturbative series can now be recursively calculated. Expected values of any (classical) observable 
f{x,p,t) is given by 



(/) = 



E 



a„ (x, t) J {x, y, t)* y" exp 



dxdy 



(24) 



where f{x,y,t) = J f {x,p) e~'^^ . This way all predictions supported on the classical reahsm can be obtained 
although satisfying quantum dynamics exactly as present on the first two coefficients of the perturbative series. Higher 
order terms are also influenced by quantum dynamics but it is not clear that in the full phase space reconstruction 
interferences cancel out in order to recover the classical character of the theory. If this was the case averages calculated 
by eqn(p^ will result on the values expected by classical theories. However it is noticeable that this doesn't happens 
and concerning eqn(^^ some results from Quantum Theory are equivalent to the present formulation of the Stochastic 
Mechanics. For example the expected value of the (classical) Hamiltonian is, using eqns(pj) and(18) 



a2 {x, t) dx + 



2toA2 



Oq (x, t) dx + / V (x, t) gq {x, t) dx 



The value of a2 is obtained by direct integration of eqn(15) and use of the field equation for 5*, the Shroedinger 
equation. The result is 



02 {x, t) ^ - 



dx^ 



9x2 



dx 



Qq ix,t) 
2A2 



(25) 



Using the condition that norm of the probability amplitude is equal to one the integral of eqn ( |25| ) is 

-h^ f , . , -n^ f _a2* . ft2 



02 (x, t) dx 



2m 



92* 
^* dx 

dx^ 



2toA2 



resulting that (H) = (H) j^jq exactly, where (H) j^jq is the expected value of the Hamiltonian calculated by Quantum 
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Mechanics for the same problem. In particular for stationary states of bound systems the classical energy becomes 
quantized. Consequently in the context presented here it appears that Stochastic Mechanics is a larger class of 
dynamical problem than Quantum Mechanics which seems be a subset of the former. In fact Stochastic Mechanics 
gives a maximum possible information in phase space through a true independence between position and momentum 
variables that is not present in Quantum Theory. 

But the above results show also that we have a disagreement between both theories. Remember that Wigner 
function in phase space satisfies Liouville equation only approximately whereas our probability density satisfy 
it exactly for zero noise power. More generally Wigner function for pure states in the adjoint space is given by 
W (x, y,t) = ^> (^x ~ i?/, t) 4' (x + iy, i)* and admits a series expansion of the form 



Wix,y,t) = ^ayjn (x,t): 



The first two coefficients of this series coincide with ours given by eqns (17) and (^ but the value of the second order 
coefficient given by 



aw2 (x, = ^ 



*(x,t) 



d^^ix,t)* 

(9x2 



+ ** (x,<) 



92^(0;, ^) 

9x2 



(9* (x,t) 



dx 
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differ from the corresponding coefficient 02 {x,t) calculated in this work, even for infinite vacuum wavelength, and 
shown in eqn (|25|). This conclusively shows that Wigner function W {x,p,t) and $ {x,p,t) are not equal. 

Having the terms in the perturbative series next step is phase space reconstruction. This task is not trivial because 
involves the whole series above displayed but if momentum of the mechanical system under study is much larger than 
Pvac = fi/^ the gaussian weight may induce a strong convergence in the perturbative series and truncation has some 
sense. Thus we expect rapid convergence for sufficiently high energy of the mechanical system compared to the mean 
vacuum fluctuation energy £. In this case the result up to second order is 



The above result shows that as an approximation, negative phase space probability also exist in our theory as it 
happens with Wigner function in Quantum Mechanics. Of course exact probability distribution should be non- 
negative and depends on calculation of all terms in the series for x- It may happens that covergence be not uniform 
a fact that complicate experimental verification since an oscilating series presents intrinsic measurements dificulties. 
Interestingly more, a new effect rises out from eqn(26) as shown in its exponential prefactor where a low energy 
phenomena controU a high energy cutoff in phase space, being this a solution somewhat different from those usually 
found to eliminate ultraviolet divergences in field theory. Even though of fundamental importance this effect will not 
considered here for further discussions. 



In this work a non-dissipative Brownian motion was considered and a perturbative series for the associated Fokker- 
Planck equation presented. Noise source was assumed to be vacuum fiuctuations of some quantum field that couples 
to particle's classical dynamics. Surprisingly the first two terms of this series accommodates quantum dynamics 
quite naturally with properties and structure of its Hilbert space stuff. Nevertheless our approach is purely classical, 
explicit in the physical meaning of all mathematical objects treated here. Moreover our probability density function 
and Wigner function differ suggesting that both theories can be compared at a experimental level specially if phase 
space tomography measurements with massive particles becomes effectiveness. 

In one line of reasoning this may indicate that our approach do not describe a truly quantal phenomena since 
a disagreement with Wigner function occur. But an alternative way may be given. Wigner function was originally 
constructed as a mathematical object without any explicit classical equivalence in what concern its (classically-defined) 
measurement. Its simultaneous dependence on position and momentum variables indicates it lives naturally in phase 
space. In addition it also satisfies approximately Liouville equation but since it can present negative values no support 
exists for an interpretation it is the quantum version of phase space probability density function. 

Quite differently, in the formalism presented here phase space dynamics is real by construction and the interpretation 
of all mathematical objects are philosophically self-consistent since they satisfy all paradigms of Classical Mechanics 
reasoning. Even the probability waves and the associated wavelengths y have clear physical meaning. In addition 
phase space reconstruction shown in eqn(26) presents in its lower order expansion negative probability values thus 
disclaiming the interpretation where this phenomenon is a mark of the Quantum Theory. In fact it is known from 
long that small noise perturbation series of the FPE may present negative terms Q. These considerations demands 
that experiment should definitively test whether dynamics in phase space as presented in this work is correct or not. 

The most serious objection against the formulation of Stochastic Mechanics we present is probably the possibility of 
a material body gain some energy from vacuum. This reasoning breaks the naive classical paradigms we claim to give 
to the rest of the theory. It is the price paid for get so smoothly quantum dynamics from a purely classical theory. We 
have argued that although vacuum power seems unphysical at first sight may in fact be a real phenomenon with true 
detection chance. If real is reasonable think that dynamcs of elementary particles are then correctly described by an 
appropriate FPE as discussed here. Vacuum temperature must be so small that usual quantum phenomena behave 
like high temperature ones and this in our interpretation means that typical wavelengths of the probability density 
are much smaller than vacuum wavelength so the whole perturbation series can be approximated by its stationary 
Liouville solution. Consequently important quantum results like atomic spectra can be understood inside our classical 




(26) 



C. Concluding Remarks 
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reasoning. Moreover our formalism open tips to think about experimental detection of alternative theories that not 
only embraces Quantum Mechanics as well test its foundations at a more precise level. These considerations are, in 
our assessment, sufficient to take account in serious way the possibility that non-dissipative Brownian motion driven 
by vacuum fluctuations be a real phenomena turning out Stochastic Mechanics a faithful description of Nature. 
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